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DISPOSITIF ET PROCEDE AUTOMATIQUE DE CALIBRATION 
GEOMETRIQUE D'UN SYSTEMS D 1 IMAGERIE PAR RAYONS X 

La presente invention concerne un dispositif et tin 
procedS automatique de calibration g&ometrique d'un 
systSme d 1 imager ie & rayons X. Des systemes de ce type 
sont utilises principalement dans le domaine medical 
5 dans le but d 1 analyser des structures internes de corps 
humains notamment de structures vasculaires. Le but de 
I 1 invention est de permettre, par une calibration aisee 
du systSme, son utilisation pour Svaluer 
quantitativement, et non plus uniquement 

10 qualitativement, les dimensions des structures internes 
analysees. L» invention se rapporte essentiellement aux 
systemes dont le detecteur de rayonnement est 
bidimensionnel, par exemple, sensiblement plan. Elle se 
rapporte, pour le procSde, aux systemes tournants 

15 d' imager ie par rayons X. 

D ans l e c^_jilun,_SYst6me d' i mager ie B. rayons X 

compose d'une source de rayons X et d*un detecteur 
bidimensionnel , tournant bu non f 1' operation geometrique 
intervenant dans la production de l 1 image est une 

20 projection conique d , un objet analyst, deploye dans un 
espace S trois dimensions (3D) , sur un espace Si deux 
dimensions (2D) qui est celui du plan de projection 
correspondant au plan de detection. Les parametres 
geometriques deer ivant completement la pro j ection 

25 conique sont au nombre de neuf. Ces neuf parametres 
peuvent etre presentes selon differentes 

parametrisations . 

Une des parametrisations possibles est la suivante: 
- trois premiers parametres concernent les trois 

30 coordonnees de la source s de rayons X dans le repSre 
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R de I'espace 3D dans lequel est plac6 I'objet 
analyst ; 

— trois autres paramfetres concernent Les trois angles 
d'Euler du changement de repere, entre le repere R de 
I'espace 3D de I'objet et un repdre R 1 li£ Sl la 
projection conique f c'est a dire forme d»un repfere 2D 
du plan de projection P et d'un troisieme axe 
orthogonal au plan de projection; 

- trois derniers parametres concernent les trois 
coordonnSes d f un point arbitraire I du plan de 
projection P (par exemple le centre de 1" image) dans 
le repere R* 

Pour simplifier les calculs dans 1' invention on 
verra par la suite qu'il y a lieu de preference de 
retenir des coordonnees cylindriques pour les trois 
premieres et trois dernieres coordonnees. Le repere R 
est alors un repere cylindrique. 

La connaissance de tout ou partie de ces parametres 
est tres souvent utile en radiologic, notamment 

lorsqu^une donnee quantitative su r un _ob jet 3D_ est 

estimee a partir d'une mesure faite dans la projection 
2D, ou dans plusieurs projections 2D, de cet objet, Des 
param&tres de chaque point de vue doivent alors etre 
connus. Un point de vue concerne 1 f orientation du 
systSme par rapport a I'objet. A titre d' exemple , on 
retiendra le probleme classique de l 1 estimation de la 
section d»un vaisseau &{ partir d f une ou plusieurs 
projections 2D, Or il est souvent impossible, ou trop 
imprecis, d'acceder a ces paramStres de fa?on directe, 
c'est a dire, par exemple, en mesurant directement sur 
le systSme d 1 acquisition la distance entre la source de 
rayons X et le detecteur. 

On appelle calibration geometrique d'un systeme 
d 1 imager ie, l 1 operation qui aboutit a la connaissance 
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indirecte precise des parametres geometriques qui 
interviennent dans la production d'une image. Le 
principe, classique, est base sur !• utilisation d'un 
fantome geometrique connu dans I'espace 3D, et dont on 
5 acquiert la projection 2D. La suite des operations 
effectuees dans ce but comporte les etapes suivantes: 

- on dispose d'un objet connu, un fantome de 
calibration pre sent ant un certain nombre de points 
caracteristiques dont la position dans 1'espace est 

10 connue par des coordonnees mesurees par rapport Sl tin 

repSre propre a cet objet; 

- on acquiert 1* image de ce fantQme dans les 
conditions geometriques d'un point de vue (ou 
incidence) que l'on veut calibrer; 

15 - on reconnalt les projections des points 

caracteristiques dans l 1 image. Pour cela on associe 
chaque point caracteristique de 1* objet a sa trace 
dans l 1 image acquise, la projection; 

- on inverse au sens mathematique le systeme 
20 d» equations decrivant la projection; 

- et on obtient finalement l v ensemble des paramStres 
de la projection pour le point de vue donne. 

Dans l f etat de l'art, la reconnaissance des points 
caracteristiques dans la projection du fant6me de 

25 calibration, est faite "manuellement" par un operateur 
humain. Par exemple, une forme de fantome de calibration 
geometrique souvent utilisee est celle d'un cube, aux 
buit coins duquel sont disposees des billes metal liques 
opaques aux rayons X. Parfois des billes supplementaires 

30 sont rajoutees pour augmenter la precision geometrique 
de la calibration. Or, la projection au sens des rayons 
X produit une image "en transparence" de 1* objet, dans 
laquelle, selon 1 1 or ientation du cube, il peut etre trds 
difficile d'associer sans se tromper les traces 2D des 
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billes avec leurs homologues 3D. On peut se ref&rer, 
pour l'Stat de I'art de la calibration g€ometrique en 
radiologie aux articles suivants: 

(1) D.L. Parker, J.Wu, D.L Pope, R. Vein Bree f G.R. 
5 Caputo and H.W. Marshall, "Three-dimensional 

reconstruction and flow measurements of coronary 
arteries using multi view digital angiography" , in New 
Developments in Quantitative Coronary Arteriography, 
J.C. Reiber and P.W. Serruys Eds., pp. 225-247, Kluwer 
10 Academic Publishers, 1988; 

(2) D.J. Hawkes, A.C.F. Colchester and C.R. Mol; "The 
accurate 3-D reconstruction of the geometric 
conf iguration of the vascular trees from X-ray 
recordings," in Physics and Engineering of Medical 

15 Imaging, R. Guzzardi Ed., Nijhoff, 1987. 

(3) M. Garreau, J T L Coatrieux, R. Collorec and 
Chardenon, "A knowledge-based approach for 3-D 
reconstruction and labeling of vascular networks from 
biplane angiographic projections", IEEE Medical 

_20 Imagingv-vol^llV-tt^-T--^ 

Cette necessite d 1 intervention humaine est un 
inconvenient majeur, fet peut m§me devenir prohibitive 
dans certains cas, lorsque le nombre de points de vue & 
calibrer est important. C'est le cas en particulier pour 

25 tout systSme d 1 acquisition de projections 

bidimensionnelles par rayons X dans le but d'une 
reconstruction 3D. Un tel systeme de reconstruction est 
par exemple decrit dans la demande de brevet f ran?ais 
N° 2 644 590 dSposSe le 20 mars 1989. Quelle que soit la 

30 m€thode de reconstruction 3D, il est necessaire au 
pr€alable de connaitre parfaitement les paramdtres 
geometriques qui caracterisent chacune des projections • 
Dans le cas probable oH une mSthode de calibration 
indirecte est necessaire, il semble indispensable 
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d'effectuer cette calibration de fa?on automatique. 

La pr§sente invention propose done un dispositif 
permettant une calibration automatique , qui s 1 applique a 
tout systeme d 1 imager ie bidimensionnelle a rayons X, et 
5 en particulier a tout syst&ae utilisS a des fins de 
reconstruction 3D* 

Le dispositif de 1* invention est tel que la 
reconnaissance, sans Equivoque, des correspondances 
entre les billes et leurs traces dams les projections 

10 est automatique. Ceci est obtehu en choisissant un 
fantSme dans lequel les billes sont distributes, de 
proche en proche, en une succession telle que des 
altitudes de billes, mesurtes le long de I'axe de 
rotation du systeme d 1 imager ie, et surtout un axe du 

15 fantome, soient monotones croissantes (ou decroissantes) 
avec un numSro d'ordre des billes dans la succession, 

Dans ce but l f invention a done pour objet un 
systSme d^magerie par rayons X avec un dStecteur 
bidimensionnel situe en vis a vis d'un tube a rayons X, 

20 ce dStecteur- et- ce— tube pouvant dans ^^ertains— cas 

tourner ensemble autour d f un axe, ce syst&me d 1 imager ie 
par rayons X comportant tin dispositif de calibration 
geometrique muni d'un fantSme de dimensions connues et 
de moyens de mesure pour deduire, a partir d f images en 

25 projection du fantdme sur le detecteur bidimensionnel, 
des coefficients de calibration de ce systSme d'imagerie 
pair rayons X, caracterise en ce que le fant6me comporte 
une succession de proche en proche de structures 
cellulaires a absorption radiologique dif f trente, plus 

30 elevSe ou plus faible, par rapport a leur environnement, 
chaque structure cellulaire etant identifiable 
automatiquement par une caracteristique ordonnee 
d 1 identification de cette structure cellulaire. 

Un autre probleme a rtsoudre pour l 1 invention est 
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liS au nombre des calibrations & entr epr endr e . Un 
systdme d 1 imager ie tournant de ce type est en effet 
susceptible d'occuper 1000 positions angulaires, ou 
plus, sur un tour. Des calibrations pour toutes les 
5 positions angulaires, meme menies automat iquement, sont 
trop longues a effectuer pour pouvoir Stre rSpetfees 
pSriodiquement, toutes les semaines ou tous les mois. II 
faut en fait prSvoir line journee de travail. 

Le procede de 1* invention vise a r^soudre cet autre 

10 problSme en faisant l^ypothtdse que, si les dimensions 
theoriques de construction du syst&me d f imagerie ne 
peuvent pas Stre considSrSes comme respectees dans le 
syst&me reel, on peut neanmoins admettre que le systeme 
d'imagerie est 4 rotation circulaire parfaite ou presque 

15 parfaite. L'idSe du proc€d§ est alors de rechercher les 
parametres intrinseques du systeme: en definitive ceux 
qui donnent en plus la direction de son axe de rotation. 
A partir de ces parametres intrinseques , qu'on peut 
connaitre trSs prScis&nent, on peut ensuite calculer 

20 pour - chaque — position angulaire.. des coefficients.- de 
calibration. On montre que, par exemple, une centaine de 
calibrations conduisent 4 une connaissance trSs precise 
de ces parametres intrins&ques. Avec cette connaissance, 
on calcule 900 autres "calibrations". On en connait 

25 alors 1000 au total en n'en ayant mesurees que 100. Le 
gain de temps est considerable. 

Ii 1 invention a done egalement pour objet un procede 
de calibration geometr ique d • un systeme d 1 imager ie par 
rayons X a d^tecteur bidimensionnel tournant autour d f un 

30 axe, en vis a vis d'un tube a rayon X et compoirtant les 
etapes suivantes: 

— on place un fantome de dimensions connues entre le 
tube et le detecteur, 

- on mesure, pour une position «en rotation de 
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1" ensemble tube dStecteur par rapport au fant6me, dans 
l 1 image radiologique du fantdme projetSe sur le 
dStecteur bidimensionnel, des coordonnees de lieux 
images de points caracteristiques du fantdme, 

- on en dSduit, pour cette position en rotation de 
1* ensemble tube dStecteur, et dans un repere associe & 
la position du fantome, des coefficients de calibration 
relatifs aux positions respectives d'un foyer de 
rayonnement du tube a rayons X et du d^tecteur 
bidimens ionnel r 

- et on r&it&re ces deux derniSres etapes pour des 
positions desirees differentes de 1' ensemble tube 
detecteur, 

caracteris4 en ce que 

- on traite les coefficients de calibration relatifs 3. 
un certain nombre de positions en rotation de 
1" ensemble tube detecteur pour en extraire des 
parametres intrinseques du systeme d f imagerie par 
rayons X qui sont indSpendants de la position en 
rotation de l f ensemble tube detecteur, „ 

- et on deduit de ces parametres intrinseques des 
coefficients de calibration relatifs & une position 
quelconque en rotation de 1» ensemble tube detecteur en 
evaluant les valeurs de fonctions analytiques de ces 
param&tres intrinseques pour une valeur angulaire de 
position en rotation de l 1 ensemble tube detecteur. 

L 1 invention sera mieux comprise £ la lecture de la 
description qui suit et £ l'examen des figures qui 
1 1 accompagnent . Celles-ci ne sont donnees qu'S titre 
indicatif et nullement limitatif de l f invention. Les 
figures montrent: 

- Figure 1: un systeme d*imagerie utilisable pour 
mettre en oeuvre 1» invention avec la presentation des 
differents reperes d f acquis it ion; 
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- Figures 2a a 2d: respectivement des vues de profil 
et de face du fantdme selon 1 1 invention, la trace des 
billes dans les projections, et un detail de mise en 
place des billes par incrustation; 

5 - Figure 3: la mise en place du fantdme dans le 

systdme d 1 imager ie; 

- Figures 4a et 4b: des images en projection d'une 
forme particuliere du fantdme apr&s traitement; 

- Figures 5a a 5d: des representations geometriques 
10 permettant une meilleure comprehension du procede; 

La figure 1 montre un systeme d'imagerie utilisable 
pour mettre en oeuvre l 1 invention. Ce syst&me comporte 
un tube 1 a rayons X avec un foyer S de rayonnement. Un 
detecteur 2 de rayonnement, par exemple du type 

15 amplif icateur de luminance utilise en radiologie, est 
place en vis a vis du tube 1. Une normale au plan P de 
detection du detecteur 2 passe sensiblement par le 
foyer S et le centre I de ce detecteur. Ceci peut §tre 
une necessite pour la formation d 1 images classiques, 

20 mais n'est. pas directement lis au probl&me^ de_ la 

calibration. En particulier, on peut calibrer, avec 
l 1 invention, un tomographe du type a translation 
homothetique. L' ensemble peut toumer avec un m§canisme 
(non reprSsente) de type connu pour les 

25 tomodensitom&tres, ou avec un arceau, autour d'un axe 3 
montre ici verticalement mais qui, dans la pratique, est 
normalement horizontal. Dans les cas des tomographes, il 
n'y a que des translations et un centre d»homothetie, 
sans meme n^cessairement un axe. II en est de meme bien 

30 stir des appareils statiques. Le systdme comporte 
egalement un processeur de traitement 4 qui fait subir 
aux signaux electriques mesures par le detecteur 2 un 
traitement contenu dans un programme 5 stocke dans vine 
memoire programme 6. Des images, ou des rSsultats 
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statistiques resultant de ces traitements, peuvent Stre 
visualises sur un moniteur 7. 

L 1 invention permet de rSsoudre le probleme de la 
reconnaissance automat ique des points homologues 2D/3D, 
5 en utilisant comme f antome de calibration un f antome 
particulidrement adapte montre sur les figures 2a a 2d 
et 3. 

Ce f antome est dans un exemple, un cylindre 8 de 
plexiglass, qui sert de support a une serie de petites 

10 billes metalliques incrustees dans la parol du cylindre , 
et disposees rSgulierement selon une helice. Les billes 
ont toutes un diametre identique, par exemple 1 mm, sauf 
la bille du milieu , qui a un diamStre plus Sieve que les 
autres, par exemple 3 mm. Les param&tres de l'hSlice 

15 sont gross ierement adaptes au systeme d' acquisition des 
images, de fa9on k ce que sa projection soit la plus 
grande possible sur le dStecteur tout en restant 
contenue dans l 1 image. De mdme, le pas de 1» helice et le 
pas d'echantillonnage des billes sur I'hSlice sont 

20 choisis de telle sorte qu'un nombre suffisant de- billes 
apparaisse dans la projection, sans qu'elles puissent se 
superposer. Les figures la a Id sont cotees en 
millimetres • Pour eviter toute deformation du fantdme 8, 
notamment par vieillissement du plexiglass le cylindre 8 

25 est obturS des deux cotes par un couvercle 9 qui en 
assure la rigiditS. 

Le fantome 8 est adapte par exemple 3. un syst&me 
d 1 acquisition const itue d f une source de rayons X et d'un 
detecteur de champ rond de diamdtre 37 cm, situ6 a une 

30 distance d f environ 135 cm, alors que l f objet est a 
environ 100 cm du foyer S. Ces chiffres sont des 
exemples qui ne sont en rien restrictifs, de meme que 
les dimensions du fantSme indiquees en annexe. Dans 
l 1 exemple il y a 61 billes. Une des coordonnSes des 
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billes, !• altitude mesur§e sur un axe longitudinal z, 
6volue de proche en proche de 5 mm en 5 mm avec une 
precision de 5/100dme. II y a 20 billes reparties 
regulierement par tour de l'h£lice. 
5 Les trois particularity principales du fant6me 8 

sont les suivantes. D'abord le fantome 8 dispose d'un 
nombre suffisant de billes pour la calibration. Elles 
sont bien reparties dans l'espace 3D. Leurs projections 
sont egalement bien reparties dans le plan P (condition 
10 "sine qua non pour une calibration 3D). La rSgularite de 
la distribution n'est cependant pas indispensable. Dans 
ce but les billes peuvent se numeroter f aci lenient en 
fonction de leur cote selon l'axe z. Ceci est obtenu ici 
selon une formule trds simple: 
15 z = h.n+zo, 

z = coordonnees selon Z, 
zO = coordonnee de la bille numero 0, 
h - pas d'Schantillonnage des billes selon Z, 
n = numero de la bille. 
20 C'est en fait la relation d'ordre qui nous importe ici: 

si zl>z2, alors nl>n2 
zl, z2 — coordonnSes des billes 1 et 2 selon Z 

nl, n2 = numero des billes* 
Si le fant6me est positionne de fa?on judicieuse 
25 devant le detecteur de rayons X, c'est a dire 
essentiellement de telle fagon que l'axe Z soit a peu 
pres parallele au plan P du d&tecteur 2, et si la 
conicite du faisceau n f est pas trop importante, cette 
relation d'ordre se retrouve egalement dans les 
3 0 projections 2-D. 

En effet, figure 3, soit V l'axe dans le plan P de 
projection qui correspond Sl la projection de l'axe Z. On 
aura la relation suivante: 

si vl>v2, alors nl>n2 
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vl, v2 = coordonnees des projections 
des billes 1 et2 selon V 
nl f n2 = numeros des billes 
Dans la suite , et par souci de simplicity , on 
5 supposera que l'axe V coincide sensiblement avec 1'axe 
vertical dans les images de projection, c'est a dire est 
parallele aux colonnes de l r image, si ce n'est pas le 
cas, on peut introduire une rotation dans le plan de 
projection pour s»y ramener. 
10 Toujours par souci de simplicity pour la suite, on 

dSfinit le repSre dans l'espace 3D comme etant life a 
l^elice: c f est a dire que 1'axe Z de l^felice dSfinit 
un des axes du repdre R, les deux autres axes Stant 
perpendiculaires a Z. Dans ce repere, les coordonnees 
15 des billes s'expriment tres facilement. Ceci signifie 
que, une fois la calibration faite, le syst&me possdde 
un axe vertical qui est celui q^occupait l'axe du 
fant6me au cours de la calibration. Des reconstructions 
3D relatives a des calibrations differentes ne sont done 

2 0 pas directement comparables . Ceci sera par contre 

possible avec le procede de l f invention vu plus loin. On 
a represents dans la figure 3 le fantome 8 positionne 
entre la source S de rayons X, et le detect eur D de 
rayons X. Le fantome est positionne de telle fagon que 
25 son axe principal Z est a peu pres parallele au plan P f 
et tel qu'il se projette selon l'axe V du detecteur. 

Enfin les billes peuvent se numeroter de fagon 
absolue. MSme si le haut et le bas de la projection de 
l'helice sont tronquees (c»est a dire que l»on ne voit 

3 0 pas la premiSre bille ni la dernifere bille du fantome) , 

on a une reference absolue qui est la grosse bille 
particulidre du milieu. La bille particuli&re, ou deux 
billes particuliSres , ou autres, sont done de preference 
placees dans le fantome a un endroit qui a de grandes 
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chances d'etre dans le champ: au milieu du fantome. 

Pour qu'on comprenne le principe mis en oeuvre dans 
le procSdd de calibration, on a deer it dans les 
paragraphes precedent un cas particulier de fantdme qui 
rSpond aux critSres exposes. En fait, tout fantSme ayant 
les particularities suivantes permet une calibration 
automat ique. Ces particular ites sont: 

- Le fantdme doit disposer de structures cellulaires, 
pour simplifier des billes, visibles en rayons X, et 
bien rSparties dans l'espace 3D. Leur contraste 
radiologique par rapport §. leur environnement doit 
etre £lev£. Les projections des billes doivent etre 
contenues dans 1 1 image 2D au moins pour les billes 
centrales. II n'est pas du tout necessaires que le 
support des billes soit un cylindre creux de 
plexiglass. 

- II doit exister une relation d'ordre entre un 
numero affectable aux billes et leurs coordonnees, 
aussi bien dans l'espace 3D selon I'axe Z, que dans la 
pro j ection selon un axe, par jexemple JL • axe V, et 
quelle que soit 1 f orientation du fantdme dans 
l'appareil. L'axe du fantame doit etre sensiblement 
parallSle au plan du detecteur. La solution des billes 
r€parties reguliSrement sur une helice comme dans le 
cas particulier du fantome decrit, n'est pas du tout 
la seule. 

- On doit pouvoir reconnaitre au moins une des billes 
parmi les autres, de fa9on & permettre une 
numerotation absolue meme si on ne voit pas la 
premiere et la derniere bille dans la projection. 
Cette bille de reference peut se distinguer par autre 
chose que par son diametre, par exemple sa forme, son 
opacity aux rayons X, son absence, ou encore une 
irregular ite du pas de distribution, ou autre. 



2700909 



13 

A titre d f exemple, on citera une autre realisation 
possible pour le fantome: les billes sont reparties sur 
vine helice "degener€e" , qui est portee par une sphere de 
plexiglass au lieu d'etre portee par un cylindre. Dans 
5 toute la suite , on raisonne avec le fantome pr§f6r6 
"helice" tel que decrit ci-dessus. Toute la mfethode 
exposee ci-apres s 1 applique egalement avec un autre 
fantome repondant aux critSres. 

On indique ci-apr&s la suite des operations a 
10 effectuer pour calibrer le systeme. 

1. Le fantdme 8 decrit ci-dessus est positionnS de 
telle fa?on que son axe principal (axe Z) soit Sl peu 
pres paralldle au plan du d&tecteur. C'est le cas par 
exemple lorsqu'on le positionne couche sur une table 
15 d'examen d , un syst&me angiographique, et que I'arceau 

qui porte le tube a rayons X et le detecteur, est 
positionnS dans un plan perpendiculaire a la table. 
Dans cette position, les projections des billes 
verifient la relation d'ordre ci-dessus. 
20 . _ 2. Une projection numerisee du fantdme est acquis e 
dans les conditions d 1 incidence que l'on veut 
calibrer. 

3. L 1 image est transf&r£e au processeur 4 de 
traitement d 1 image. 

25 4. L 1 image est corrig^e des distorsions geometriques 

Sventuelles, par exemple avec une methode telle que 
celle dScrite dans la demande de brevet f rangais 
N° 2 633 793 deposee le ler juillet 1988. 
5. L* image corrigee I est ensuite soumise a un 

3 0 algorithme de detection, localisation, et 

labellisation automatique des billes decrit ci-aprSs. 
6* Les coordonnees des billes dans le rep&re 3D, R, et 
de leurs projections dans le repere 2D du plan P, sont 
exploitees pour estimer les param&tres geometriques de 
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la projection: les coefficients de calibration. Uh 
algorithme base sur une minimisation d'erxeur par la 
technique du gradient conjugue est propose et 
egalement decrit ci-aprSs dans ce but. 

On decrit ici un exemple possible de !• algorithme 
de calibration, qui permet de passer de l 1 image en 
projection du fantSme de calibration, figure 
4a - image I, £ des paramStres geometriques intrinsSques 
de la machine. On montre tout au long de cette 
description les images resultant des differents 
traitements . 

1. D f abord on segmente 1* image I. L'objectif est ici 
de supprimer dans 1 9 image I le signal de la projection 
du cylindre de plexiglass qui supporte les billes. On 
utilise pour cela des outils de morphologie 
mathematique: transformation chapeau haut de forme, 
avec un element structurant de forme ronde et de 
taille superieure a la taille des billes dans 
l 1 image I. On obtient une image SI qui ne contient 
plus que le signal correspondent a l\absprptiqn_ des 
billes. Les outils de morphologie mathematique sont de 
type connu et publies dans de nombreux ouvrages. 

2. L ■ image SI doit etre purifiee des artefacts 
residuels, provenant des bords du cylindre. On utilise 
pour cela de nouveau une transformation de morphologie 
mathematique visant a supprimer les structures 
verticales dans 1» image. L 1 image resultat est une 
image S2 visible figure 4b en negatif. L 1 image S2 
contient done toutes les billes (signal positif ) sur 
un fond a valeur nulle, au bruit pres. 

3. On segmente ensuite l 1 image S2. I^objectif est ici 
de ne garder que le signal provenant de la grosse 
bille centrale du fant6me. On utilise pour cela une 
erosion de l 1 image S2 # d*une taille superieure a la 
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taille des petites billes mais inf€rieure a la taille 
de la grosse bille. On obtient vine image S3, qui ne 
contient que la grosse bille* 

4. On d€tecte ensuite toutes les billes dans l 1 image 
5 S2. On veut ici localiser au pixel prSs les 

projections de toutes les billes dans I 1 image. On 
commence & la ligne 0 de 1" image (en haut du 
d£tecteur) . On parcourt toute l 1 image du haut vers le 
bas. Un signal positif rencontrS dans 1' image est 

10 declare etre une bille si ce signal est superieur Sl un 

certain seuil 1. On donne le numero 0 a la premi&re 
bille rencontree, puis le numSro est increments 3. 
chaque bille rencontree, jusqujen bas de 1» image. On 
est sGr que I'ordre des numeros est le bon, grSce aux 

15 propriety de l'hSlice dScrite ci-dessus. 

5. On calcule ensuite les coordonnSes exactes des 

m 

billes. L'objectif ici est de connaltre les 
coordonnees de chacune des billes dans 1* image, et 
ceci avec une precision sub-pixel. Pour chaque bille 

20 dans . 1» image S2,- on .effectue un calcul de centre^ de 

gravite en niveaux de gris de I 1 ensemble des pixels 
autour de la bille detectee, et tels que leur niveau 
de gris soit superieur a un deuxieme seuil, seuil T2 
(seuil T2 < seuil Tl) . Un calcul de ce type est par 

25 exemple entrepris dans la deuxieme demande de brevet 

frangais citee ci-dessus. Son principe repose sur le 
fait que les billes sont suffisamment fines pour ne 
pas absorber tout le rayonnement X. Alors, en 
l 1 absence de saturation d f absorption une precision 

30 inferieure au pixel d' image est obtenue. 

6. On detecte ensuite la grosse bille dans 1» image S3; 
on connait alors ses coordonnees au pixel pres. On 
compare ses coordonnees avec eel les de toutes les 
billes dStectSes dans l 1 image S2, et on reconnalt 
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ainsi dans la liste de toutes le bille celle qui est 
la grosse bille- 

7. La grosse bille etant reconnue dans la projection, 
on peut maintenant renumSroter toutes les billes de 
5 fa?on absolue, et ainsi Stre assure d f une 

correspondance exact e, pour chaque bille, entre son 
numero et ses coordonnees 3D, connues par ailleurs, et 
les coordonnees 2D de sa projection* 

Connaissant, pour chaque bille, ses coordonnSes 3D 
10 et les coordonnees 2D de sa projection, on peut calculer 
1» erreur de reprojection associSe & un choix de valeurs 
arbitraires pour les 9 param&tres de la projection 
conique. Cette erreur est obtenue par la difference 
entre la position calculee avec ces valeurs arbitraires 
15 et la position trouv&e dans 1* image du fantame de 
calibration. On peut done estimer les paramStres ou 
coefficients de calibration associ£s au point de vue, 4 
la position en rotation, que I'on veut calibrer en 
choisissant les parametres minimisant un critere 
20 d f erreur de reprojection pour l 1 ensemble des billes 
detectees. 

L 1 algorithme propose d 1 estimation des parametres 
utilise la technique de minimisation du gradient 
conjuguS applique au critere d 1 erreur quadrat ique 

25 moyenne. Cette technique classique necessite le calcul 
du critSre et de ses dSrivees premieres en f onction des 
param&tres. Cet algorithme est public aux noms de W.H. 
Press, B.P. Flannery, S.A. Teukolsky and W.T. 
Vetterling, " Numerical Recipes in C", Cambridge 

30 University Press, 1988, chapitre 10.6, page 317 et 
suivantes. D'autres algorithmes, decrits au mSme endroit 
dans cet ouvrage, et denommes Simplex et Newton sont 
aussi utilisables. Le premier, Simplex, 4 l'avantage de 
s» adapter a n'importe quel systSme coherent de 
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param&trisation. II ne converge cependant pas 
rapidement. L'algorithme du gradient conjugue converge 
plus vite mais necessite un choix de parametrisation tel 
que les expressions analytiques des paramStres soient 
5 d€rivables. C'est le cas de la param§trisation prSsentSe 
ci-dessus. Les calculs relatifs a cet algorithme on ete 
presences au congrds "3D Advanced Imaging Processing in 
Medicine" a Rennes en France en 1992 au cours de la 
14eme conference annuel le IEEE EMBS, du 2 au 4 Novembre 
10 1992 , par Anne Rougee, Catherine Picard, Cyril Ponchut 
et Yves Trousset. 

Les etapes en sont les suivantes: 

1. Initialisation: les paramStres doivent etre 
initialises a des valeurs realistes, ce qui permet de 

15 lancer la minimisation dans une region proche de la 

solution. 

2. Normalisation: autour de cette valeur initiale on 
peut calculer la valeur d'un crit&re correspondant 3. 
un deplacement d'une units sur chague axe de l'espace 

20- des paramdtres, puis effectuer un changement. de 
variable dans cet espace par changement d'echelle sur 
chaque axe afin que tout d§placement unitaire sur 
chaque axe provoque le meme effet sur le critere. 

3. Minimisation: partant de la valeur initiale, et 
25 tenant compte d'un changement d'dchelle trouve au 

cours de l*etape de normalisation, on applique la 
technique iterative du gradient conjuguS pour 
minimiser 1 1 erreur quadrat ique moyenne sur l 1 ensemble 
des billes detectees. A chaque iteration on calcule 
30 une nouvelle estimation des parametres en fonction des 

deriv€s du crit&re. Le processus est stoppe lorsque 
l^cart relatif entre 1* estimation courante et la 
suivante est infSrieure a un seuil donn§. 

On appelle machine de reconstruction 3D une machine 
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qui effectue des radiographies d'un objet sous 
diff erentes incidences obtenues en tournant autour d'un 
axe. Ces radiographies sont ensuite utilises pour 
reconstruire une information 3D. Cette machine peut par 
5 exemple etre construite a base d'un anneau rigide sur 
lequel tourne !• ensemble tube detecteur auquel cas l'axe 
de rotation est fixe et ne change jamais au cours du 
temps. Ou bien elle peut etre construite a base d'un 
arceau inclinable supportant I 1 ensemble tube detecteur , 

10 auquel cas l'axe de rotation peut etre positionne de 
diff erentes fagons dans 1'espace. 

Plut6t que de calculer pour chaque point de vue f 
pour chaque position en rotation de 1' ensemble tube 
detecteur, les param&tres dits coefficients de 

15 calibration, on peut selon un proc§d6 de 1' invention, 
calculer des parametres intrinsSques du t syst£me 
tournant. Ces parametres intrinsdques sont alors des 
fonctions analytiques de la position en rotation de 
1' ensemble tournant. II suffit ensuite, pour obtenir les 

20 parametres de calibration relatifs a une position^ de 
calculer la valeur des parametres intrins&ques pour une 
valeur donnSe d 1 or ientation du systSme tournant. Les 
systdmes de l'art connu delivrent deja normalement une 
mesure, le plus souvent sous forme d'un signal 

25 electrique, repr^sentant leur orientation dans 1'espace. 
II suffit d'utiliser cette mesure pour obtenir les 
coefficients de calibration relatifs a cette orientation 
en la portant dans les expressions analytiques des 
param&tres intrinsdques. 

30 Pour estimer les parametres intrins&jues decrivant 

sous forme condensee la gSometrie globale d'une machine 
de reconstruction 3D, on peut proceder de la fa?on 
suivante: estimer les parametres de la projection 
conique, vue par vue, lors d'une premiSre etape, puis 
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estimer les parametres intrins&ques & partir de 
l 1 ensemble de ces parametres lors d'une deuxieme Stape. 
Premiere etape 

Dans le cas d'une machine de reconstruction 3D, la 
5 methode de calibration automatique est particuliSrement 
bien adapt&e pour estimer les coefficients de 
calibration vue par vue au cours d'une rotation. 

On positionne l»h€lice sur la table d'examen, son 
axe Z approximativement centre sur l'axe de rotation de 
10 la machine. Pour chacun des points de vue d 1 acquisition 
de la machine au cours d f une rotation, on va acquerir la 
projection de l'hglice, et traiter ainsi toutes les 
images acquises, avec la methode dScrite prec^demment. 

La paramStrisation proposee - avec des coordonnSes 
15 3D cylindriques - est adaptee a l 1 estimation des 
parametres le long d'une trajectoire circulaire autour 
de I'axe Z. I*es &tapes d' initialisation et de 
normalisation sont effectuSes pour la premiere vue 
seulement. Les valeurs des paramdtres trouvSs pour la 
2 0 premiere vue sont ensuite utilisfees pour- initialiser 
l f estimation des parametres de la vue suivante et ainsi 
de suite jusqu'S. la dernier e vue. Ceci permet 
d'accelerer les calculs. 

On obtient done les paramStres geometriques pour 
25 chaque point de vue, qui sont stockes sur le disque 
numerique de la machine, et serviront dans toute etape 
de reconstruction 3D. 
Deuxieme etape 

II est possible selon l 1 invention de passer de la 
30 connaissance des coefficients de calibration, vue par 
vue, de la machine, a la connaissance de parametres 
intrinseques dScrivant sous forme condensee la gSom&trie 
globale de la machine. 

Soit N le nombre de points de vues calibres. Par 
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exemple N vaut 100* Dans le cas oH chaque point de vue 
est calibre independamment de son voisin, on dispose de 
9N parametres pour decrire la geomStrie des 
acquisitions. Or dans le cas d , un systeme tube detecteur 
5 lies,, qui effectue une rotation autour d»un axe fixe, on 
voit bien que ces 9N parametres sont redondants. La 
geometrie d' acquisition de la machine de reconstruction 
3D peut se parametrer avec un nombre beaucoup plus 
limite de parametres, qui sont; de preference: 
10 - les coordonnees d f un isocentre O du systSme 
(3 parametres) ; 

- la direction de l'axe Z» de rotation passant par 
l 1 isocentre O. Elle est donnee par 2 angles 
(2 parametres) ; 

15 - les rayons rj et rg des trajectoires de I et de la 

source S (2 parametres) ; 

- La difference de cote (suivant l'axe de rotation 
Z 1 ) entre la trajectoire de la source et celle de I 
(1 parametre) ; 

20 - l^cart angulaire e entre 1* angle polaire de_S par 

rapport 5 Z» et 1' angle polaire de I par rapport a Z 1 
(1 parametre) ; 

- les 2 tilts ai et a 2 du detecteur, les angles 
d* orientation du ddtecteur par rapport a l'axe Z 1 

25 (2 parametres) ; 

- 1 1 angle c*3 entre les colonnes de 1 1 image et la 
projection de l'axe Z 1 dans I 1 image (1 parametre) 

On obtient done un total de 12 paramdtres. D*autres 
choix pour ces douze parametres sont possibles. Celui-ci 
30 est prefSrS pour des considerations analogues aux 
precSdentes: simplification des calculs. On notera que 
dans ce cas le refSrentiel dans lequel sont exprimees 
les valeurs de ces parametres intrinseques est 
independant de la position du fantome. II est done 
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tr an spo sable d'une calibration a l 1 autre. Une methode 
pour estimer ces 12 parametres en fonction des 9N 
parametres de depart peut etre par exemple la suivante: 
1* Trouver la direction de I'axe de rotation Z f : 

On considere les N points I n (n=l,N) formant la 
trajectoire de I. En thSorie, cette trajectoire est 
plane et orthogonale a Z 1 • On cherche done le plan le 
plus proche, au sens des moindres carr§s, de l 1 ensemble 
des I n et on en d&duit la direction de Z 1 , qui est la 
normale a ce plan. 

Soient (X n ,Y n , Z n ) les coordonnees d f un point I n . 
On cherche le plan d 1 Equation z = ax + by + c qui soit 
le "plus proche" de 1* ensemble des I n . On appelle 
d(a,b,c) la distance entre le plan z = ax + by + c et 
l f ensemble des points In, On recherche les coefficients 
(a,b,c) qui minimisent la distance d(a,b,c) • On 
definit: 

d(a,b,c) = 2 (ax n + by n + c - z n ) 2 
et on trouve alors tres facilement le triplet (a,b,c) 
qui minimise d(a,b,c) par regression linSaire. 

2, Trouver les coordonnees de l'isocentre, et le rayon 
de la trajectoire de I: 

L'isocentre est le centre a la fois de la 
trajectoire de I et de celle de S. On va done rechercher 
ce centre en analysant la trajectoire de I. 

Remarque: dans cette etape, on cherche les 
parametres (centre et rayon) de la trajectoire de I, 
puis dans l f etape suivante on cherchera le rayon de la 
trajectoire de S. On aurait pu, en toute logique, faire 
1» inverse , c f est & dire commencer par chercher le centre 
et le rayon de la trajectoire de S, puis en dSduire le 
rayon de la trajectoire de I. Cette solution est 
toutefois mo ins robuste que celle exposSe ici, car le 
point S etant plus eloigns que le point I de I'axe de 
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rotation, les informations sur la position de S sont 
plus entachees d'erreur que celles sur la position de I. 
On utilise done les informations "stlres" (celles sur X) 
comme base de la methode et non pas celles sur S f qui 
sont plus incertaines. Une autre solution consisterait 
egalement a trouver de manidre simultan^e les rayons et 
centres des trajectoires de I et de S. Pour les memes 
raisons de robustesse, cette methode n'a pas StS 
retenue . 

Comment trouver les coordonnees de l^socentre, 
e'est-a-dire du centre de la trajectoire de I? On 
dSfinit un plan n comme etant orthogonal a Z' et ayant 
une cdte suivant Z» Sgale a la moyenne des cdtes des 
points I n . On projette orthogonalement sur ce plan, en 
des points I n tous les points I n . La trajectoire de I 
etant supposee circulaire, les l' n sont a priori situes 
sur un cercle. On cherche alors le cercle le plus 
proche, au sens des moindres carres, de l 1 ensemble des 
I'n ce qui nous fournit: 

- le rayon ri de la trajectoire de I (e'est le - rayon 
du cercle que 1 1 on a trouv6) f 

- les coordonnees de I'isocentre, qui est ici defini 
comme 1* intersection de rr et de l'axe de rotation 
(e'est le centre du cercle que l'on a trouvfe) • 

Comment definir et trouver le cercle le plus proche? 
Soit C(A,r) un cercle de centre A et de rayon r. Ce 
qu*il faut, c»est dSfinir tine distance d(A,r) entre le 
cercle C(A,r) et ^ensemble des I'n* Plusieurs mSthodes 
sont possibles. 
Methode 1 - figure 5a 

Soient (x'nfY'n) les coordonnees de I'n* ° n 
d€finit: d(A,r) = S I'n^'n le point J f n etant defini 
comme 1 9 intersection du cercle C(A,r) et de la droite 
Al'n- 



2700909 



Methode 2 - figure 5b 

Quand la machine de reconstruction 3D emet des tirs 
X suivant des positions angulaires reguli&res, c'est §L 
dire tous les SB degres (59 constant) , on peut utiliser 
5 cette connaissance pour dtre plus robuste dans la 
determination des parametres du cercle. 

On va dSfinir d(A,r) par: d(A,r) = Z I' n J f n , le 
point J 9 n etant defini cette fois-ci comme le point du 
cercle C(A,r) qui possede, par rapport & une direction 
10 de reference, un angle polaire: e = 9q + (n-1) * 58. 

Que I'on adopte la methode 1 ou la methode 2 pour 
def inir la distance entre le cercle et le nuage de 
points, la suite est la meme : on recherche, par une 
methode quelconque de minimisation (gradient conjugue 
15 par exemple) les parametres A et r du cercle (ou/ pour 
la methode 2, les paramStres a, r et Bq) qui minimi sent 
la distance au cercle. Une fois trouve ces paramfetres, 
c'est fini, puisque : 

- Le rayon de la trajectoire de J" est connu : c'est r 
20 - L'axe Z' de rotation est entierement connu, puisque 
l'on a determine son orientation a l'etape precedent et 
que l'on sait maintenant que cet axe passe par A. 
3» Trouver le rayon de la trajectoire de S 
On connait le centre de la trajectoire de S, 
25 puisque c'est par definition le meme que celui de la 
trajectoire de I, que l'on vient de trouver i. l'Stape 
precedente. Cherchons maintenant le rayon • On considere 
les points S' n , projections des points S sur le plan 7T. 
Le rayon de la trajectoire de S est dSfini comme la 
30 moyenne des distances entre l f isocentre O et les points • 
S' n : 

rayon de I = l/N* Z0S' n 

4. Trouver la difference de cSte suivant I 1 axe Z' 
entre la trajectoire de la source et celle de I. 
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C'est la moyenne des cotes suivant l'axe Z* des 
points S n : voir figure 5c. On en deduit la difference de 
cote par la distance au plan 7r. 
5. Trouver I 1 angle e: 
5 On travaille a partir des projections I* et S' n 

dans le plan n. Pour chaque n, on calcule les angles 
polaires (par rapport a O) de I' n et S'n et on appelle 
e n leur difference (voir figure 5d) . On estime e come 
etant la moyenne des e n . 

10 6. Trouver les tilts ai et oe2 du detecteur 

Pour chaque vue, on calcule dans le repere lis & 
l'axe de 1'hSlice, les coordonnees de la normale au 
detecteur, Ensuite, par application d'une formule de 
changement de repSre, entre le repere lie & l'axe de 

15 l'h§lice et celui li§ h 1'isocentre et a l'axe de 
rotation Z • , on en deduit les coordonnees de la normale 
au detecteur dans ie repere li& 1 1 isocentre et a 1 * axe 
Z' de rotation. Cette formule de changement de repdre 
est deduite de la connaissance du rep&re lie a l'h&lice. 

20 Par ail leur s on connait 1' isocentre puis que c'est le 
centre du cercle trouvS pr ecedemment . On connait 
egalement 1 ' orientation de 1 ' axe Z 1 . Les deux autres 
axes de ce repSre sont figes 3. une rotation pr6s. De 
preference on fait en sorte, pour les f iger, que la 

25 premiere vue soit affectee dans le nouveau repdre d'un 
angle Bq lie a une position privilegiee de la machine. 
Dans ce nouveau repere, ces coordonnees ne dependent que 
des 2 angles c*i et 0:2- Pour chaque n f on calcule ces 2 
angles notes a± n et 0^2x1* Ensuite ai est defini comme la 

3 0 moyenne des axh et idem pour c*2n* 
7. Trouver 1' angle a.3 

Pour chaque vue, on calcule dans le repere lie a 
l'axe de l'helice, les coordonnees du vecteur unitaire 
parallele aux colonnes du detecteur. Ensuite, par 
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application de la formule de changement de repere, on en 
dSduit ses coordonn€es dans le repere lie a l'isocentre 
et t l*axe Z 1 de rotation, or, dans ce nouveau rep&re, 
ces coordonnees ne dependent que de l 1 angle 0*3. Pour 
chaque n on calcule cet angle, note <*3 n . Ensuite 0*3 est 
defini comme la moyenne des 0^3 n . 
Troisieme etape 

On connait maintenant 12 paramdtres intrinsdques 
qui definissent la geometrie du systeme. A partir de la 
connaissance de ces 12 parametres intr insSques , il est 
trfes facile de recalculer, pour chaque vue n, les 9 
paramStres locaux de calibration d^finissant cette vue. 
Ce calcul est d f autant plus simple qu'il a lieu dans le 
repSre intrinseque de la machine, c'est-a-dire dans le 
repere lie a l^socentre et a l'axe de rotation. On 
appelle (rs n , as n , zs n , theta n , phi n , psi n , ri n , ai n , 
zi n ) les 9 parametres locaux relatifsa la vue n, tels 
qu'ils sont d^finis dans Particle publie au Congres de 
Rennes. II vient alors : 

rs n _ = rayon de la trajectoire _de la source 

(= const ante) 

as n = e 0 + 180 + (n-l) * SB 

zs n = difference de cote calculee prScSdemment 
(= constante) 

theta n = 9 0 + (n-l) * SO + al 
phi n = a2 
psi n = a3 

ri n = rayon de la trajectoire de S (= constante) 

ai n = 9q + (n-l) * SQ 

zi n = 0 (= constante) 

Ces neuf parametres ou coefficients de calibration 
sont done determines a partir de 1 ■ incidence d§duite de 
e 0 + (n— 1) £9. lis sont connus avec beaucoup de 
precision compte tenu des evaluations statist iques dont 
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lis resultent. Avec des mesures precises de 1» incidence 
on arrive alors S. des coefficients de calibration 
ca leu les tres precis. L 1 incidence elle meme est d£duite 
de l , indice n. Cet indice est connu par la division de 
la circonf erence de la machine en un nombre donne de 
positions prSdtablies. 

On notera que, de preference, on recalculera, a 
partir des expressions analytiques selon les 
coefficients intrinsdques, les coefficients de 
calibration relatifs aux N vues calibrees, De ce fait, 
les coefficients de calibration pour ces N (100) vues 
sont evalues avec autant de precision que pour les 
autres vues (les 900 autres) . Le resultat d 1 ensemble est 
encore meilleur. 

Une autre mSthode possible consisterait Si estimer 
directement les paramStres globaux de la machine & 
partir de l 1 ensemble des coordonnees des billes 
dStect£es dans toutes les projections de l»helice. Mais 
cette methode est longue. 
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REVENDICATIONS 



1 - Syst&me d' imager ie par rayons X avec un 
detecteur bidimensionnel situe en vis §l vis d'un tube k 
rayons X, ce systeme d f imager ie par rayons X comportant 
un dispositif de calibration geometrique muni d'un 
5 fantome de dimensions connues et de moyens de mesure 
pour d^duire, a partir d 1 images en projection du fantdme 
sur le detecteur bidimensionnel , des coefficients de 
calibration de ce systSme d'imagerie par rayons X, 
caractSrisd en ce que le fant6me comporte une 

10 succession, de proche en proche, de structures 
cellulaires Si absorption radiologique contrastSe par 
rapport a leur environnement , chaque structure 
cellulaire Stant identifiable automatiquement par une 
caractSristique ordonnSe d 9 identification de cette 

15 structure cellulaire, 

2 - Systeme d' i mager ie par rayons X selon la 

revendication 1, caracterisS en ce que, ce dStecteur et 
ce tube tournant ensemble autour d'un axe, le fantome 
comporte Sgalement un axe, un support, et une 

2 0 distribution des structures cellulaires sur ce support, 
de telle fa?on que des altitudes de ces structures 
cellulaires, mesurSes selon cet axe du fantdme, soient 
monotones et croissantes avec un numSro d'ordre de ces 
structures cellulaires, I'axe du fantome etant destine a 

25 etre orients dans le systSme d* imager ie par rayons X, en 
phase de calibration, sensiblement paralldlement 4 I 1 axe 
de ce systdme d 1 imager ie par rayons X, et de telle fa£on 
que cette relation d'ordre se retrouve dans la 
projection du fantdme* 

30 3 - SystSme d 1 imager ie par rayons X selon la 
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revendication 2, caractSrist en ce que le support est un 
cylindre, de preference circulaire, dont les 
generatrices sont paralleles a l*axe de ce cylindre et 
en ce que la succession ordonnee des structures 
5 cellulaires sur ce support a sensiblement la forme d , une 
helice. 

4 - Systdme d'imagerie par rayons X selon la 
revendication 3, caracterise en ce que le support 
cylindrique est muni a ses extremites de plaques de 

10 calage pour empScher sa deformation. 

5 - SystSme d'imagerie par rayons X selon l*une des 
revendications 1 a 4, caractSrisS en ce que la 
caractSristique ordonnee d 1 identification de chaque 
structure cellulaire est constitute par un repSrage de 

15 cette cellule par rapport a au moins une structure 
cellulaire particuliere du fantome. 

6 - SystSme d'imagerie par rayons X selon la 
revendication 5, caracterist en ce qu'une structure 
cellulaire particuliere du fantome possSde une taille ou 

20 vine absorption — radiologique — differente de- - ce-lle- des 

autres structures cellulaires et est placSe sensiblement 
au milieu de la distribution des structures cellulaires. 

7 - Systdme d 1 imager ie par rayons X selon l'une des 
revendications 1 a 6, caracterisS en ce que les 

25 structures cellulaires possSdent des absorptions 
radiologiques insuf f isantes pour empScher totalement le 
passage du rayonnement X a travers elles. 

8 - Syst&ne d'imagerie par rayons X selon I 1 une des 
revendications 1 £ 7, caracterise en ce que les 

30 structures cellulaires comportent des billes metalliques 
placees au fond de puits realises dans un support. 

9 - Precede de calibration geometrique d»un systSme 
d , imagerie par rayons X a detecteur bidimensionnel 
tournant autour d'un axe, en vis a vis d'un tube a rayon 
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X et comportant les etapes suivantes: 

- on place un fantome de dimensions connues entre le 
tube et le detecteur, 

- on mesure, pour une position en rotation de 
I 1 ensemble tube detecteur par rapport au fantdme, dans 
1' image radiologique du fantdme projetSe sur le 
dStecteur bidimensionnel, des coordonnees de lieux 
images de points caracteristiques du fantdme, 

- on en deduit, pour cette position en rotation de 
l 1 ensemble tube detecteur , et dans un repere associe a 
la position du fantdme, des coefficients de 
calibration relatifs aux positions respectives d'un 
foyer de rayonnement du tube a rayons X et du 
detecteur bidimensionnel, 

- et on rSitSre ces deux dernieres Stapes pour des 
positions dSsirees diffSrentes de l f ensemble tube 
detecteur, 

characterise en ce que 

- on traite les coefficients de calibration relatifs a 
un certain nombre de positions en rotation de 
1' ensemble tube dStecteur pour en extraire des 
parametres intrinsdques a2, <*3...) du systSme 
d*imagerie par rayons X qui sont independants de la 
position en rotation de l 1 ensemble tube d&tecteur, 

- et on dSduit de ces parametres intrinsSques des 
coefficients de calibration (rs n , as n ...) relatifs a 
une position quelconque en rotation de l f ensemble tube 
detecteur en evaluant les valeurs de fonctions 
analytiques de ces parametres intrinsSques pour une 
valeur angulaire de position en rotation de l f ensemble 
tube detecteur. 

10 - Procede selon la revendication 9, caracterise 
en ce que pour deduire des coefficients de calibration 
relatifs a une position desiree en rotation de 
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1* ensemble tube ddtecteur, 

- a) on simule, avec des coefficients de calibration 
arbitraires, des coordonnSes de lieux theoriques qui 
soient des images , sur le detecteur bidimensionnel, de 

5 points caracteristiques du fantdme, 

- b) on compare les coordonnees de ces lieux 
theoriques avec celles des lieux mesures dans l 1 image 
du fantdme, 

- c) on modifie en correspondance les coefficients de 
10 calibration arbitraires, 

- d) et on reitdre les deux Stapes a) et b) jusqu'a ce 
que la comparaison montre un difference nSgligeable, 
auquel cas on retient comme coefficient de calibration 
pour cette position en rotation de 1 1 ensemble tube 

15 detecteur les derniers coefficients de calibration 

modifies. 

11 - Procede selon la revendication 10/ caractSrise 
en ce que 

- on remplace les coefficients de calibration 
.20 arbitraires_ par__ des coef f icients„ de calibration 

deduits de la construction mecanique du systeme 
d r imager ie par rayons X. 

12 - ProcSdS selon la revendication 10 ou la 
revendication 11, caract§ris6 en ce que 

25 - pour dSduire des coefficients de calibration 

relatifs Si une nouvelle position desiree en rotation 
de 1" ensemble tube detecteur, decalee par rapport a 
une position en rotation precedente, on utilise S la 
place des coefficients de calibration arbitraires 

30 relatifs a cette nouvelle position en rotation, les 

derniers coefficients modifies obtenus pour cette 
position en rotation precedente. 

13 - ProcSde selon l f une des revendications 9 £ 12, 
caractSrisS en ce que les coefficients de calibration 
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sont exprimes par des fonctions analytiques des 
paramdtres intrinseques , ces fonctions analytiques etant 
decrites dans un repere cylindr ique . 

14 - Proced£ selon l'une guelconque des 
revendications 9 a 13 caracterise en ce que les 
parametres intrinseques sont obtenus par evaluation 
statistiques & partir d'un nombre restreint d v images du 
fantSme pour diver ses orientation du systeme d» imager ie. 

15 - ProcSde selon I'une quelconque des 
revendications 9 a 14 caracterise en ce que on d6duit, Si 
partir des parametres intrins&ques, des coefficients de 
calibration relatif s a des positions en rotation de 
1» ensemble tube detecteur pour lesquelles on a effectu6 
des mesures des coordonnees de lieux images de points 
caracteristiques du fantdme, et pour lesquels on en 
avait dfeduit des coefficients de calibration • 
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